Anomalous roughening in competitive growth models with time-decreasing rates of 

correlated dynamics 



o 

(N 
D 

m 



o 



in 

in 
o 

O 

m 



X 



F. D. A. Aarao Reis 
Instituto de Fisica, Universidade Federal Fluminense, 
Avenida Litoranea s/n, 24210-340 Niteroi RJ, Brazil 
(Dated: February 5, 2013) 

Lattice growth models where uncorrelated random deposition competes with some aggregation 
dynamics that generates correlations are studied with rates of the correlated component decreasing 
as a power law. These models have anomalous roughening, with anomalous exponents related to 
the normal exponents of the correlated dynamics, to an exponent characterizing the aggregation 
mechanism and to that power law exponent. This is shown by a scaling approach extending the 
Family- Vicsek relation previously derived for the models with time-independent rates, thus providing 
a connection of normal and anomalous growth models. Simulation results for several models support 
those conclusions. Remarkable anomalous effects are observed even for slowly decreasing rates of 
the correlated component, which may correspond to feasible temperature changes in systems with 
activated dynamics. The scaling exponents of the correlated component can be obtained only from 
the estimates of three anomalous exponents, without knowledge of the aggregation mechanism, and 
a possible application is discussed. For some models, the corresponding Edwards- Wilkinson and 
Kardar-Parisi-Zhang equations are also discussed. 

PACS numbers: PACS numbers: 81.15Aa, 05.40.-a, 05.50.-|-q, 68.35.Ct, 68.55.-a 



I. INTRODUCTION 

The technological interest in thin films, multilayers and 
related nanostructures motivated the study of continu- 
ous and atomistic growth models and kinetic roughening 
theoriesi"— . In the simplest cases, surface fluctuations fol- 
low the standard Family- Vicsek (FV) scaling relation^i^. 
The amplitudes of fluctuation modes saturate at times 
increasing with the wavelength and the global roughness 
saturates when the largest excited mode is of the order of 
the substrate size. Local and global surface fluctuations 
scale with the same exponents. This is called normal 
roughening. On the other hand, a large number of sys- 
tems show anomalous roughening^^— , where local slopes 
continuously increase and local exponents differ from the 
global ones^'^'^. Theoretical approaches have shown var- 
ious conditions, such as symmetries and conservation 
laws, where anomalous scaling (AS) is expected^i^. 

With the advance in deposition techniques, changing 
physico-chemical conditions during the growth of a film 
is possible. In vapor methods, this can be achieved by 
varying pressure or temperature. In electrodeposition, 
the evolution of surface morphology changes the electric 
field distribution and, consequently, the local reaction 
rates may change. Thus, local growth is a competition 
between different dynamic processes (adsorption, reac- 
tion, surface diffusion) with time-dependent rates. As a 
simple illustration, one may consider a constant flux of 
atoms or molecules and time- varying rates of surface pro- 
cesses (responsible for fluctuation correlations). These 
features may have significant effects on roughness scal- 
ing. For instance, in Ref J^, it was shown that a sudden 
change in growth conditions of an Edwards- Wilkinson^® 
model may lead to power-law relaxation between different 
growth regimes. Moreover, Ref shows that stochastic 



growth equation models with time-dependent couplings 
have AS under certain conditions. 

Those findings motivate the present investigation of 
lattice models with competition of two aggregation dy- 
namics with time-dependent rates. We consider limited 
mobility models where uncorrelated random deposition 
competes with some type of correlated component whose 
rate decreases in time as a power-law. These systems 
show AS, in contrast to the normal scaling observed in 
the original models with constant rates^^. The anoma- 
lous exponents are related to the (normal) exponents of 
the correlated growth dynamics and to the exponent of 
the time-decreasing rate. Those relations are obtained by 
direct substitution of the time-dependent rate in the FV 
relation of the model with time-independent rates, which 
provides a connection between scaling relations of normal 
and anomalous models. Those conclusions are supported 
by simulation results. Remarkable anomalous effects are 
observed even for slowly decreasing rates, which may cor- 
respond e. g. to small changes of temperature in systems 
with activated dynamics. For some models, the general 
forms of the corresponding Edwards- Wilkinson (EW)^ 
and Kardar-Parisi-ZhangiS, equations are also discussed. 

The hydrodynamic equations associated to our growth 
models have time-dependent couplings, thus the exis- 
tence of AS is predicted by the approach of Reffi^. How- 
ever, the aim of work is not discussing conditions for AS, 
but to explore details of a class of lattice models with 
realistic dynamics showing that feature, with a simple 
connection of normal and anomalous roughnening and an 
important role of the aggregation mechanism. A possible 
application of our model is discussed, as well as systems 
where it is certainly not applicable. 

The rest of this work is organized as follows. In Sec. 
II, we review the scaling concepts for competitive models 
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FIG. 1: (Color online) Illustration of the model with com- 
petition between RDSR (probability p) and UD (probability 
1 — p) in one-dimensional substrates. In (a), after the choice 
of RDSR, the incident particle attaches at the column at the 
right side. In (b), after the choice of the RDSR, the particle 
randomly chooses between the neighboring columns (curved 
arrows); only the attachment at column at the right side is 
illustrated. In (a) and (b), after the choice of UD, the particle 
sticks at the top of the column of incidence. 



with correlated and uncorrelated components and extend 
the approach to the FV relation of the height-height cor- 
relation function. In Sec. Ill, we extend that approach 
to competitive models with time-decreasing rates of the 
correlated component and calculate their anomalous ex- 
ponents. Sec. IV shows simulation results of some models 
in two-dimensional substrates, confirming the theoretical 
predictions. In Sec. V, we discuss the EW equations 
(in all dimensions) and KPZ equations (in 1 -I- 1 dimen- 
sions) associated to some of those models. In Sec. VI, we 
discuss the relation with experimental work and the re- 
liability of our models. Sec. VII summarizes our results 
and conclusions. 



II. LATTICE MODELS WITH COMPETITION 
OF CORRELATED AND UNCORRELATED 
DEPOSITION 

In these models, cubic particles of size a are sequen- 
tially deposited on randomly chosen columns of a d- 
dimensional substrate of lateral size L, in a direction 
perpendicular to the substrate. The time interval for de- 
position of one layer of atoms [{L/a)'^ atoms] is r. Each 
incident particle has two possibilities: (1) permanently 
stick at the top of the column of incidence, with proba- 
bility 1—p (uncorrelated deposition - UD); (2) execute 
some movement (diffusion, desorption, or eventually stick 
at some point of the column of incidence) according to an 
aggregation rule that takes into account the neighboring 
column heights, with probability p. This second compo- 
nent generates correlations of the column heights and is 
hereafter called correlated deposition (CD). 

Fig. 1 illustrates the model where the correlated 
component is random deposition with surface relaxation 
(RDSR)i^. In RDSR, the particle sticks at the top of 



the column of incidence if no neighboring column has 
a smaller height, otherwise it sticks at the top of the 
column with the smallest height among the neighbors 
(in case of a draw, one of the columns with the small- 
est height is randomly chosen). The competitive model 
UD-RDSR and another one involving UD and ballistic 
deposition (BD)^ were introduced by Albano and co- 
workers2ii22. 

In this section, we consider the case where the proba- 
bility p is kept constant during the film growth. 

The global roughness of the film surface is defined as 
the rms fluctuation of the height h around its average 
position h: 



W{L,t)= Uh-h) 



1/2 



(1) 



where the overbars indicate spatial averages and the an- 
gular brackets indicate configurational averages. At short 
times, UD dominates, thus the roughness increases as 



Wrd ~ a{t/T) 



1/2 



(2) 



After a crossover time ic, the CD determines the uni- 
versality class of the process, characterized by a set of 
exponents hereafter indexed with C . The roughness fol- 
lows Family- Vicsek (FV) scaling^ as 



W{L,t) « AL°'^ 



(3) 



where ac is the roughness exponent, / is a scaling func- 
tion such that / ~ 1 in the regime of roughness saturation 
(t -T- oo) and tx is the characteristic time of crossover to 
saturation, which scales as 



(4) 



where zc is the dynamic exponent. For t <^ (but 
t 3> tc), the roughness scales as 



(5) 



where f5c = ctc/zc is the growth exponent. In this 
growth regime, f{x) ^ x^'^ in Eq. ([3]). 

The exponents ac, Pc and zc depend on the basic 
symmetries of the CD, but the amplitudes A, B and C 
arc modcl-dcpendcnt. For small p, they scale as 



and 



p 



B 



C^p- 



(6) 



(7) 



(8) 



where the convention of crossover exponents {5, y, 7) of 
Albano and co-workers^^ was used. FV scaling implies^ 



y/3c - S + '-f = 0. 



(9) 
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Different scaling approaches have aheady explained 
the values of the exponents obtained in simulations of lat- 
tice modelsiii^^. Below we review the scaling approach 
of Ref.^'' for the global roughness and extend it to the 
height-height correlation function (HHCF) and the local 
roughness. 

For p 1, most deposited atoms attach to the top 
of the randomly chosen column (UD). Thus, the height 
difference of neighboring columns, Ah, is of order Ah ^ 
a(Ai/r)^^^ after a time interval At (as Eq. [2]). On the 
other hand, the average time for a CD event (probability 
p) to take place at a given column is ~ r/p. 

If CD is BD, that event immediately creates correla- 
tions between the neighboring columns. Thus, the time 
of crossover from UD to correlated growth is 

^ Tc - p^V {BD - like). (10) 

This applies to other BD-like models, as discussed e. g. 
in Ref.24. 

However, if CD is a solid-on-solid (SOS) model, such as 
RDSR, a single correlated event does not cancel the ran- 
dom fluctuation of neighboring column heights Ah. In- 
stead, that fluctuation is cancelled only when the number 
of correlated events Nc is of order Ah/ a. At the crossover 
time tc, that number is Nc = tc/Tc tcp/r. Thus 

tc^p-\ {SOS) (11) 

In both cases, all time scales of the purely correlated 
system {p = 1) are also changed by the scaling factor 
tc/r, such as the saturation time tx (Eq. 2]). Thus we 
have ?/ = 1 for BD-like models and y = 2 for SOS models 
in Eq. ([7]), independently of the class of the CD and sub- 
strate dimension (this improves the picture that emerged 
from a former work^^, which restricted y = 1 to KPZ and 
?/ = 2 to EW). 

The average height difference between neighboring 

1 /2 

columns saturates at Ahc ~ a{tc/T) , which is of order 
ap-y^^ both for BD-like and SOS models. This is the 
scaling factor for global height fluctuations (Eqs. [3] and 
[S]), consequently 

6 = y/2. (12) 

It gives (5=1/2 for BD-like models and 5 1 for SOS 
models. For this reason, a single exponent {y) fully char- 
acterizes the UD-CD crossover. 

If one is interested in comparison with experimental 
data, two other quantities are more useful. The first 
one is the local roughness w (r, t) , which is averaged over 
windows of size r gliding through the surface. The second 
one, which will be the focus of most of our calculations, 
is the HHCF of columns at distance r: 

G {r, t) = {[h {rh + r, t) - h {rh,t)f),T = \r\. (13) 

Here, the configurational averages (angular brackets) are 
taken over different initial positions rp, different orien- 
tations of r and different deposits. The local roughness 
and the HHCF have the same scaling properties. 



For the competitive models analyzed here, the HHCF 
in the regime of thin film growth {tc <^tx) scales as 

The constants A and B scale as the corresponding ones 
for the global roughness (Eqs. [6] and [7]) and g is a scaling 
function. An alternative form that helps connection to 
anomalous scaling relation&iS is 

./WJ)^Dt^cG^-^y (15) 

where 

i5^p-y(i/2-/3c). (16) 

For smah distances r (r < {t/Bf'"''^), the HHCF is 
hereafter called local HHCF and scales as 

./G^c^p'y'\"^. (17) 

This quantity is time-independent in systems with nor- 
mal FV scaling, but time-dependent with anomalous 
scaling. For large distances, the HHCF is hereafter called 
saturation HHCF and scales as 

^ p-yi^i^-Pc)tPc _ (18) 

Note that jHc < 1/2 for any type of CD with normal scal- 
ing, thus the exponent —y (1/2 — /3c) is always negative. 

III. COMPETITIVE MODELS WITH 
DECREASING RATES OF CORRELATED 
COMPONENTS 

We consider competitive models whose correlated com- 
ponent has time decreasing probability as 

P=(^ + l) ^ (19) 

where A is positive. Eq. (|19|) is properly defined for 
lattice models, so that at t = all deposition events are 
correlated {p = 1). For long times (t ^ r), we have 

p^t-^. (20) 

Eq. ([20]) can be substituted in Eqs. and (HH) to 
give AS for the HHCF, as follows. 

The local HHCF scales with time and distance as 

y^_^,A/2^ac_ (21) 

This quantity has the same scaling of the local roughness 
in small box sizes. Thus, the local roughness exponent 
aioc is the same of the correlated component, ac, and 
the local slope exponent is 

K = yA/2. (22) 
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This exponent represents the degree of anomaly of the 
system (k = for normal growth). It depends not only 
on the exponent of the time-decreasing rate but also on 
the aggregation conditions (BD-like or SOS) through the 
exponent y. However, it does not depend on the universal 
exponents characterizing the class of the CD. 

In theoretical works, k is sometimes called the lo- 
cal anomalous growth exponent 13*^^. In experimen- 
tal works, it is frequently called /?/oc (local growth 
exponent)^, but some authors give a different definition 
for /Sjod^. For these reasons, we keep the notation close 
to that of Ref.J^, consistently with Eqs. and ((22)) . 

The saturation HHCF scales as 



/3 = /3c + Ay (1/2 -/3c). (23) 



This quantity has the same scaling of the local roughness 
in large box sizes. Here, the global growth exponent /3 
is larger than the growth exponent of the CD. In some 
works, that exponent is denoted as the sum of /3/oc and 
/3*9. 

The crossover between growth and saturation regimes 
of the HHCF occurs when r - (tx/B)^^^^ (Eq. [H]). 
This gives the characteristic time tx as 



1-yA 



(24) 



This shows that the anomalous dynamical exponent z is 
also larger than that of the CD, reflecting the decreasing 
rate of propagation of correlations. The expected relation 



; -I- aioc/z = P 



(25) 



is observed in all cases. 

Due to the continuously decreasing role of the corre- 
lated dynamics, the global roughness W (Eq. [T]) does 
not saturate, but continuously increase (i. e. there is not 
steady state). This can be shown by direct substitution 
of the time-dependent form of p (Eq. [T9)) in the FV rela- 
tion ([3]), with the predicted forms of the scaling function 
and of the p-dependent amplitudes. 



IV. NUMERICAL RESULTS 

Three applications illustrate the remarkable effect of 
time-decreasing p, even in cases where p is still of order 1 
after deposition of many layers. In all simulations, two- 
dimensional substrates {d = 2) are considered. 

First, we consider RDSR (also known as Family model 
- Fig. l)i^ as the correlated component. This model is 
in the EW class, thus (3^—0 and ac = 0, both corre- 
sponding to logarithmic increase of the roughness in time 
and size, and zc = 2 (diffusive dynamics)i£. Since it is 
an SOS model, we have y = 2, thus Eqs. and ([M]) 
give /3 = A and z = 2/ (1 - 2A). 

Simulations of the UD-RDSR model were performed 
with A = 1/6 and A = 1/4 in substrates of size L — 




1 2 
logio(r/a) 



FIG. 2: (Color online) HHCF for the RDSR (solid curves) and 
for the UD-RDSR model with A = 1/6 (dashed curves) . From 
bottom to top, deposition times are t/r = 100 (red), t/r = 
200 (green), t/r = 400 (blue), and t/r = 800 (magenta). 
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FIG. 3: (Color online) Local HHCF for the UD-RDSR model: 
ro = 8a for A = 1/6 and ro — 6a for A — 1/4. Dashed lines 
are least squares fits of the data with slopes k ~ 0.19 and and 
K 0.28, respectively. 



1024a up to t/r = 10^ In Fig. 2, we show the HHCF for 
A = 1/6 and for the pure Family model (p = 1, A = 0) 
at several times. The split of the curves for small box size 
r is a clear signature of anomalous scaling for A = 1/6. 

Fig. 3 shows the time evolution of the local correla- 
tion function (fixed r = ro) for two values of A: ro = 8a 
for A = 1/6 and ro = 6a for A = 1/4. The power-law 
increase is clear in both cases and the linear fits give es- 
timates of the exponent k slightly larger than the predic- 
tion of Eq. ([2D) for y — 2. However, exponents calculated 
from local roughness or correlation function in short-time 
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FIG. 4: (Color online) EfTective local slope exponents, with 
to — 200r, for the UD-RDSR model: upper symbols ro — 6a 
and A = 1/4, lower symbols for ro — 8a and A — 1/6. Dashed 
lines are least squares fits of the data. 



simulations of lattice models frequently deviate from the 
expected values of the model class^, thus the deviation 
in Fig. 3 is not unexpected. 

Very accurate estimates of k are obtained from the 
effective exponents 



. ^ ^ In [Gioc (r-Q , t) I Gioc (rp , t - tp )] 



(26) 



which tend to k as t — >■ cx), with fixed tp. We assume that 
the leading correction to scaling of Gioc (not of its square 
root) is a constant term, in analogy to the intrinsic width 
in the context of roughness scaling (see e. g. Ref.'??, and 
references therein). Using Eqs. (|2ip and (|22|) with y = 2, 
this assumption gives K{t) — > A with a correction term 
proportional to t^^'^. In Fig. 4, we plot K{t) as a function 
of using the same data of Fig. 3 with tp = 200r. 
The linear fits give asymptotic estimates k — 0.255 for 
A = 0.25 and n = 0.163 for A = 1/6 = 0.1666..., in 
excellent agreement with the theoretical prediction k = 
A. 

Now we consider BD as the correlated component. In 
BD, the incident particle attaches to the site where it 
has a nearest neighbor occupied site (at the sides or 
below it)^i^, creating a porous deposit due to lateral 
aggregation. The model is in the KPZ classic, where 
ac « 0.385, zc « 1.615, and (3c « 0.24^-^. Using 
y — 1 ioT the competitive model, Eqs. ((23)) and (|24)) give 
13 w 0.24 + 0.26A and z w 1.615/ (1 - A). 

Simulations of the UD-BD model were performed in 
two-dimensional substrates of size L = 1024a up to t/r = 
10^, again for A = 1/6 and A — 1/4. Fig. 5 shows the 
HHCF for A = 1/6 and for the pure BD model {p = 1, 
A = 0) at several times. Comparison with Fig. 2 shows 
a smaller split of the curves for small r, which is expected 



FIG. 5: (Color onhne) HHCF for BD (solid curves) and for 
the UD-BD model with A = 1/6 (dashed curves). From bot- 
tom to top, deposition times are t/r — 100 (red), t/r — 200 
(green), t/r = 400 (blue), and t/r = 800 (magenta). The 
long-dashed line has the expected slope of the KPZ class. 



for the smaller exponent y (Eq. [21]) . 

An interesting effect shown in Fig. 5 is the small 
slope of the curves in the scaling region, which gives lo- 
cal roughness exponents much smaller than the expected 
o^ioc ~ etc ~ 0.38. The apparent roughness exponent 
is even smaller in the competitive model. However, large 
discrepancies between exponents from simulations of BD- 
like models and KPZ values are frequent''^, thus the de- 
viations in Fig. 5 are not a particular feature of models 
with AS. 

A third model used to illustrate our results has the 
restricted SOS (RSOS) model^ as the correlated com- 
ponent. The condition for aggregation of a particle at 
the column of incidence is that the height difference be- 
tween all neighboring columns cannot exceed one lattice 
unit, otherwise the aggregation attempt is rejected. This 
means that the aggregation is accepted only if the height 
of the column of incidence is smaller than or equal to the 
heights of all neighbors (this rule is used in the exten- 
sion to the competitive model). The RSOS model also 
belongs to the KPZ class and usually shows small scaling 
corrections. Using y — 2 for the competitive model and 
estimates of KPZ exponents^^i"-^, Eqs. ([23]) and dH]) give 
/3 « 0.24 + 0.52A and z « 1.615/ (1 - 2A). 

Simulations of the UD-RSOS model are performed in 
two-dimensional substrates under the same conditions of 
the previous ones. Fig. 6 shows the HHCF for A = 1/6 
at times between t/r — 100 and t/r = 800. Again, AS 
is clearly shown by the split of the curves for small r, in 
contrast to the case of time-independent rates. 

The inset of Fig. 6 shows the time scaling of Gsat 
for the competitive model, which gives /3 w 0.35 for the 
longer times. This is in good agreement with the theoret- 
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of scaling exponents. Indeed, scaling arguments follow- 
ing the same lines of Refi^ lead to z ~ 2/(1 — fi) and 
global roughness exponent aoL = ^/ (1 ^ ^) {c^gl = 
a/oc + zffi^). Since zc = 2 and ac ~ for the normal 
EW equation (with coefficients not changing in time) and 
y = 2 (SOS model), comparison with Eq. (|24|) leads to 
H. ~ 2 A. Equivalently, the full surface tension coefficient 
is ~ (as predicted in Ref^ in 1 + 1 dimensions). 
Similar result is obtained for lattice models in the EW 
class in any other substrate dimension. 

On the other hand, the KPZ equation^^ corresponding 
to such competitive models can be constructed only in 
d= 1: 



7j{x,t). (28) 



FIG. 6: (Color online) HHCF for the RSOS model (solid 
curves) and for the UD-RSOS model with A = 1/6 (dashed 
curves). From bottom to top, deposition times are t/r = 100 
(red), t/r = 200 (green), t/r = 400 (blue), and t/r = 800 
(magenta). Inset: saturation HHCF as a function of time, 
with a linear fit of the last five points (solid curve) with slope 
13 0.35. 

ical prediction /3 w 0.33. Again, extrapolation of effective 
exponents (not shown) provide estimates much closer to 
the theoretical prediction. 

A strinking feature of this model is the small slope of 
the logG X logr curves for small r, caused by the large 
height fluctuation produced by UD, in comparison with 
the typically small roughness of the RSOS model. This 
feature masks the asymptotic scaling for small r in the 
model with AS. 



V. EW AND KPZ EQUATIONS 

For competitive models such as UD-RDSR, the general 
form of the associated EW equation can be derived, with 
coefficients related to p and scaling exponents related to 
A. Consider, for instance, the case in two-dimensional 
substrates {d — 2). The surface tension term of the EW 
equationi^ is expected to decrease in time, thus we pro- 
pose 

| = .„(|^) %^h + ,i,,t), (27) 

where h is the height at the position x in the d- 
dimensional substrate at time t, vq ^^nd to are constants, 
and ?7 is a Gaussian noise^'^* with zero mean and co- 
variance (?7 {x,t)ri{x',t')) = D5'^[x - x')5[t -t'). 

A complete representation of the lattice model would 
require calculation of the constants and i as well as 
the coefficients of higher-order derivatives not shown in 
Eq. ([?7)l . However, that form is sufficient for calculation 



The roughness amplitudes of Eqs. ^ and (jH) scale with 
the coefficients of the linear and nonlinear terms as yl ^ 
and B ~ |A|i^-i/233^ Using Eqs. dH), 0, and 
([T9| . for BD-like models, we have v t^^ and A ^ 
t"3A/2^ However, for SOS models {y = 2) in the KPZ 
class, such as UD-RSOS, we have v ^ t"'^^ and A ~ 
t"^'^. Again, a complete representation of the lattice 
model by a stochastic equation would require calculation 
of coefficients of the terms in Eq. ([28]) (not only their 
scaling properties) and of higher-order terms not shown 
there. 



VI. DISCUSSION 

The results for A = 1/6 show that anomalous scaling 
can be observed in competitive growth models even with 
very slow changes of the rates and far from the crossover 
region p <^1. ki t = 800r (Figs. 2, 4 and 5), we have 
p ~ 0.33, i. e. there is a significant contribution of the 
correlated process. For comparison, numerical work on 
competitive models with fixed p usually focus the regime 
p < 0.3i2i2ii2S. Thus, the time-dependence of the com- 
petitive dynamics is the main ingredient for the AS, in- 
dependently of approaching the crossover region. 

Although the competitive models presented here are 
drastic simplifications of a real thin film growth process, 
their assumptions are not far from experimentally fea- 
sible conditions. If one assumes that the correlated ag- 
gregation is an activated process, then the decrease of 
the probability p may be a consequence of decreasing 
the temperature during the film growth. In the case 
of an EW or KPZ model, that assumption suggests an 
Ahrrenius form for the surface tension coefficient as '^-^ 
exp (— S/ZcbT), which gives p ~ j/^/^ ~ exp{—E/ykBT) 
(y = 1 or 2). If other processes produce correlations (e. 
g. surface diffusion), then the coefficients of the corre- 
sponding growth equation will also have the Ahrrenius 
form. Now, assuming that E — 0.2eV, a temperature 
decrease from 330K to 3QQK gives a ratio of probabili- 
ties w 0.7 for y — 2. For comparison, in Figs. 2, 4 and 
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5 (A = 1/6), the probability p decreases by a factor 0.71 
from t = lOOr to t ^ 800r.' 

Experimental measurement of anomalous exponents k, 
/3 and z allows testing the hypothesis of a competition of 
uncorrelated and correlated deposition. From Eqs. p2|). 
([23|) . and (j24|) . their values provide estimates of the scal- 
ing exponents of the correlated component as 

A = 2K/y , Pc = Y^ ' zc = z{1-2k). 

(29) 

The estimates of fic and zc can be compared with ex- 
isting theories of interface growth. Note that they are 
measured independently of the exponent y, thus it is not 
necessary a priori knowledge of the aggregation mecha- 
nism. 

Applications of this model are not expected when k is 
very close to 1/2. For instance, \i n = 0.4, we obtain 
z = 5zc and /? — k = ^/Sc, i- e. the anomalous exponents 
are much larger than those of theories of normal scaling 
(the difference /? — k is frequently referred as /3 in ex- 
perimental works). For this reason, the AS with large k 
observed in metal electrodeposition by some authors has 
to be explained by other approachesM. Other systems 
showing AS with large values of ^^iiSiii, as well as exper- 
iments on metal dissolutio n^"^'^^ , will also give unreliable 
estimates of and/or zc if our competitive model is 
applied. 

On the other hand, the work on Cu electrodeposition 
by Lafouresse et ali^ suggests an application. At the 
early stages of growth and for low potential values, the 
exponents k — 0.07, /3 — k — 0.21 and a/oc ^ 1 are 
obtained. Using our model, we find /3c' = 0.24 and z = 
4.1. These values are very close to the exponents fie = 
0.25 and = 4 of the linear fourth-order stochastic 
growth equation 

^ ^iyiV^h + r]{x,t) (30) 

(Mullins- Herring class?'^^), which represents growth 
dominated by surface diffusion. Consequently, a possible 
interpretation of the experiment under those conditions 
is that the role of the surface diffusion is continuously 
decreasing in time due to a competition with an uncor- 
related random growth. 



The results for lattice models shown in Figs. 4 and 
5 also suggest that deviations in the small r scaling are 
frequent. Consequently, when experimental data is an- 
alyzed and compared to our model, one should consider 
that a disagreement in the values of a may not repre- 
sent a failure of the model. Instead, this is an expected 
feature at short times / small thicknesses—. 

VII. CONCLUSION 

We studied lattice models with competition of corre- 
lated and uncorrelated mechanisms for aggregation of 
deposited particles and where the rate of the corre- 
lated component decrease in time as a power-law. The 
Family- Vicsek scaling relation derived for the case of 
constant rates is extended by direct substitution of the 
time-dependent rate, giving anomalous scaling relations 
that are supported by simulation data. The anoma- 
lous exponents are related to the exponents of the cor- 
related growth dynamics and to the exponent of the 
time-decreasing rate. The EW and KPZ equations cor- 
responding to some models are discussed. 

Even with a slow decay of the correlated component, 
which may correspond to small temperature changes in 
systems with activated dynamics, a remarkable anomaly 
is found. Indeed, application of the model is expected 
only for systems with small anomaly exponent k (/3* or 
l^ioc in some works). This excludes most electrodeposi- 
tion works, but a case of thin Cu films electrodeposited 
at low potential values shows a possible application of 
growth dominated by surface diffusion. Since our scal- 
ing relations provide exponents of the correlated compo- 
nent depending only on anomalous exponents that can 
be measured experimentally (independenly of the expo- 
nents of the time-decreasing rate and of the aggregation 
mechanism), the test in other real systems is simple and 
may eventually help to understand their basic physical 
mechanisms. 
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